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ABSTRACT 

We present a two dimensional map of the gas temperature distribution in the Centaurus cluster, based on ASCA 
observations derived using a novel approach to account for the energy dependent point spread function. Along 
with a cool region, centered on NGC4696, asymmetric temperature variations of moderate amplitude are detected. 
The hottest region roughly coincides with the position of the second brightest galaxy NGC4709, known to be the 
dominant galaxy of one of the subgroups in the Centaurus cluster. ROSAT images show faint surface brightness 
emission also centered on this galaxy. The imaging and spectral results suggest that a subcluster, centered on 
NGC4709 is mergin g with the main cluster centered on NGC4696, in agreement with the earlier suggestion by 
Lucey et al. (1986) that these two systems are located at the same distance despite their different line-of-sight 



velocities. 

Subject headings: galaxies: clusters: individual (Centaurus) - 



X-rays: general — large-scale structure of 



1. INTRODUCTION 

The Centaurus cluster is a nearby (mean velocity of the dom- 
inant component is about 3000 km s -1 , e.g. Lucey, Currie, & 
Dickens 1986) X-ray bright cluster with a 2-10 keV luminosity 
of 7 x 10 43 ergs s" 1 (Jones & Forman 1978). The cluster X-ray 
emission is strongly peaked on the cD galaxy NGC4696, indi- 
cating a cooling flow with a moderate mass deposition rate of 
- 30-50 Mq (e.g. Matilsky, Jones & Forman 1985; Allen & 
Fabian 1994; White, Jones & Forman 1997). 

The cluster was intensively studied in X-rays with imag- 
ing instruments by Einstein (Matilsky, Jones & Forman 1985), 
ROSAT (Allen & Fabian 1994), and ASCA (Fukazawa et al. 
1994). Rather smooth X-ray isophotes, having a slightly el- 
liptical shape (Allen & Fabian 1994), and the presence of a 
cooling flow suggest that the Centaurus cluster is a relatively 
relaxed system. However, measurements of line-of-sight veloc- 
ities indicate the presence of two distinct components (Lucey et 
al. 1986), one centered on NGC4696 and another dominated by 
NGC4709, but having a factor of 2.5 fewer galaxies. These two 
components, denoted as Cen 30 and Cen 45, have mean veloc- 
ities and dispersions of v = 3041 km s" 1 , a = 586 km s" 1 and 
v = 4570 km s , a = 280 km s" 1 respectively (Lucey et al. 
1986). Although this may be the simple pr ojection of two spa- 
tially separated groups, Lucey et al. (1986) argued, on the basis 
of comparison of luminosity functions, color-magnitude rela- 
tions and galaxian-radius distributions, that both components 
are in fact located at the same distance. This would imply 
that the main cluster, Cen 30, is in the process of accreting the 
smaller group, Cen 45, at the present time. We show below that 
the ASCA observations provide strong support that the Centau- 
rus cluster is undergoing a merger. 

The structure of the paper is as follows: Section 2 briefly de- 
scribes the X-ray observations and data reduction procedure. 
In Section 3 we discuss results of the analysis and present argu- 
ments in favor of an ongoing merger of the Cen 30 and Cen 45 



subclusters. Section 4 summarizes the results. 



2. DATA ANALYSIS 

For the analysis presented below, we used three observations 
of the Centaurus cluster taken with ASCA on 1993 June 30, 
July 5 and 1995 July 19. These three observations have two 
slightly different pointing directions and together cover the re- 
gion extending from the SE to the NW. Primarily GIS data 
were used for the analysis, which has a larger field of view 
and more stable characteristics over 1993-1995, compared to 
the ASCA SIS. Standard screening has been applied to the 
ASCA data with recommended values, forming a cleaned data 
set which was used for all the analysis. Similar screening cri- 
teria have been applied to publicly available observations of 
"blank fields", which were used for background subtraction. 
The effective area was calculated using the most recent version 
of the efficiency curves (xrt_ea_v2_0.fits). A contour map of 
the resulting background subtracted, vignetting corrected, GIS 
image in the 0.5-2 keV energy band is shown in Fig. [l] (Plate 
1). The circle coinciding with the peak of the surface brightness 
corresponds to the dominant galaxy NGC4696. 

To reconstruct the projected (two dimensional) gas tempera- 
ture distribu tion, we used the approach described in C hurazov 
et al. (1996). This method provides an efficient characteriza- 
tion of the shape of the spectrum at a given position as a lin- 
ear combination of template spectra, corresponding to optically 
thin plasma emission models at two different temperatures. For 
our analysis, we used the MEKAL model from XSPEC (con- 
volved with the efficiency of the mirror systems and the ap- 
propriate ASCA detector) at temperatures of 2 and 8 keV, as 
template spectra. For both template spectra, the heavy element 
abundance was set to 0.5 of solar, the absorption column density 
to N H = 8 x 10 20 cm" 2 , and the redshift to 0.0104. The value of 
the temperature at a given position is calculated from a simple 
combination of the best fit weights of the two template mod- 
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els to the spectrum observed at this position. Due to the linear 
nature of the method no actual fitting is required and, hence, 
the method is computationally fast. Although the limitations 
of this method are significant (e.g., simple spectral forms must 
be assumed a priori), this approach does provide a convenient 
method to examine two dimensional spatial variations in gas 
temperature. Note that for a single temperature spectrum, the 
accuracy of the method is sufficient (f or temperatures above ~ 
2 keV) as shown by the simulations in |Churazov et al. (1996) 

The extended, energy dependent wings of the ASCA mir- 
rors are known to cause significant systematic errors in spectral 
analyses of extended sources, usually leading to spurious tem- 
perature increases at large radii from the cluster center where 
the scattered emission from the cluster center contributes sig- 
nificantly to the total emission. We developed an approach to 
correct for the effects of the PSF (for a detailed description in- 
cluding a more rigorous derivation see Gilfanov et al. 1998) 
which relies on the existence of a compact core in the ASCA 
PSF. For the PSF model, we use publicly available GIS images 
of Cygnus X-l observed with different positional offsets. The 
method is based on the assumption that the ASCA PSF (P) can 
be represented as the sum of two components: 

P — Pcore Pwings ( 1 ) 

where P core corresponds to the rather compact "core" of the 
PSF and P W i ngs corresponds to the extended (large spatial scale) 
wings of the PSF, in such a way that Y^Pwings "C 1, i.e., the 
wings of the PSF contain only a small fraction of the photons 
(we assume here that by definition J^P = !)■ The convolution 
of the true sky image (S) with the PSF can be rewritten as 

D = S<Z>P core + S<g>P wings (2) 

where D is the observed detector image. Provided that P core is 
sufficiently compact, the first term in the above sum, S <E> Pcore, 
would provide a resonable approximation to the sky image S, 
since it is the convolution of the true sky image with a relatively 
compact function. From (||) it follows that 

S®P core =D-S®P wings (3) 

Utilizing the fact that Y,Pwings *C 1, one can obtain an approx- 
imate solution, Y for S ® P CO re, by substituting the measured 
detector image D for S in the right side of equation (Q) 

Y = D - D <8 P wings = S <8 P core - (D - S) <g> P wings 

= S £§> Pcore ~(S(E) Pcore ~ S) (8> Pwings ~ S <8> Pwings ® Pwings (4) 

One can see that the two last terms in the above equation can 
be dropped, assuming that the cluster image S is "smooth" and 
E Pwings < 1. Thus, 

Y^S®P core (5) 

To practically implement this method, we use the following 
simple choice for P cor e and Pwings'- 

P = P- P = if r < R = 6' 

1 core 1 ■> 1 wings u ; 11 ' lK core u 
Pcore — 0) Pwings ~ P-i if ^> Pcore ~ 6 

where r is the distance from the center of the PSF. The accuracy 
of the approximation in eq. (||) is determined by the fraction of 



photons contained in the PSF wings, ; which in turn de- 

pends on the value of R cor e - the larger is R CO re, the better the 
approximation. On the other hand, for larger values of R cor e, 
the PSF core defined above becomes less compact and its en- 
ergy and position dependence becomes more important. The 
choice of R CO re = 6' was found to be adequate for clusters with 
core radii larger than a few arcminutes. 

Equation (g) involves a two-dimensional convolution (D <8> 
Pwings) of the observed detector image (in a given energy band) 
with Pwings, which in turn depends on position and energy. 
Straightforward convolution is a very time consuming proce- 
dure, but a dramatic acceleration can be achieved via a Monte- 
Carlo approach. For each detected photon, we choose N = 
10-20 random positions jc,,y,- (distributed according to P W ings) 
and subtract their contribution (J^Pwings/N) from the corre- 
sponding positions. The function P W i ngs itself is represented as a 
one-dimensional cumulative array, i.e. the two-dimensional ar- 
ray P W i ngs (of size M x M) is represented as a one-dimensional 
array of size M 2 , whose k-th element is equal to the sum of 
all elements from Pwings with indices less than k. Such a rep- 
resentation of P w tn gs allows one to generate random positions 
(distributed like Pwings) with just a few arithmetical operations, 
using the bisection method. 

The simplicity of this method makes its implementation com- 
putationally fast and it can be performed in each of the 1024 
ASCA GIS energy channels (or 512 SIS channels) on a fine 
(« 15") angular grid. On the other hand, to decrease the sys- 
tematic uncertainty intrinsic to this method, the value of R CO re 
must be sufficiently large. This restricts somewhat the applica- 
bility of this approach so that it can be applied most accurately 
to large, generally nearby clusters. In particular, the method 
has proven to work fairly well for clusters of galaxies with core 
radii exceeding several arcmin (see Gilfanov et al. 1998). 

The resulting temperature map for the Centaurus cluster is 
shown in Fig. |] (Plate 1). The statistical noise increases to- 
wards the edges of the map and the spectral measurements are 
shown only for those regions where the ratio of the upper limit 
on the temperature to the lower limit is less than a factor of 1.5. 
One can clearly see a cool region centered on NGC4696 as well 
as asymmetric variations of temperature. However, we wish to 
emphasize that these variations are of moderate amplitude with 
the temperature ranging from ^3-3.5 keV to ~ 4-4.5 keV (ex- 
cluding the cooling flow region). Note that the hottest region 
to the SE is roughly coincident with the position of NGC4709 
which is the dominant galaxy of the Cen 45 subcluster (Lucey 
etal. 1986). 

The same PSF correction method can be applied to spectra 
accumulated from arbitrary regions. We have chosen a set of 
regions as shown in Fig. || and have accumulated GIS spectra 
for these regions. Shown in Fig. || are the results of fits to the 
GIS spectra for each region with a single temperature MEKAL 
model from XSPEC (with fixed redshift of z= 0.0104, fixed in- 
terstellar absorption of = 8 x 10 20 cm" 2 ). The temperature 
of the gas, the heavy element abundance, and the normaliza- 
tion were free parameters. Note that although the actual shape 
of the spectrum may be more complex than the model we are 
considering (in particular in the cooling flow region, see e.g. 
Fukazawa et al. 1994), for our purpose this simple model pro- 
vides a useful characterization of the spectral hardness (tem- 
perature). Since we are interested primarily in the temperature 
variations outside cooling flow regions, a single temperature ap- 
proximation is sufficient. This was further confirmed by values 
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Fig. l.— Greyscale image of the gas temperature distribution in the Centaurus cluster (see also color Plate 1). The image was 
smoothed with a 2' Gaussian. Contours correspond to the 0.5-9 keV GIS surface brightness image corrected for background and 
vignetting. The surface brightness image was smoothed with a 0.5' Gaussian. Circles show the positions of the brightest galaxies 
(brighter than m=14). 
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Fig. 2. — Left: The regions selected for spectral fitting. Solid symbols indicate the bright (m < 14) galaxies in the field. Right: The 
temperature and surface brightness versus region number. Errors indicate regions of 68 per cent confidence for a single parameter of 
interest. 



of Xr ~ 1 f° r me outer regions. 

In Fig. |2|, dotted lines separate the results obtained for differ- 
ent annuli. NGC4709 falls into region #24 which has a higher 
temperature than adjacent regions. Inspection of other annuli 
shows that in general, the temperature is higher in the direction 
from NGC4696 towards NGC4709. Comparing Fig. |with the 
two dimensional temperature map (Fig. [l]) one can see that both 
methods give consistent results. 

3. EVIDENCE FOR MERGING 



As is discussed by Allen & Fabian (1994), the soft X-ray 
emission from the Centaurus cluster is slightly elliptical. They 
also noted that the second brightest galaxy NGC4709 is lo- 
cated along the semimajor axis of the cluster. We co-added 
three ROSAT PSPC images, RP800607A01, RP800192N00, 
and RP800607N00, with total exposure time of about 18000 
seconds. No cleaning was applied to the data to achieve maxi- 
mal statistical significance. Exposure maps were generated for 
each pointing and also coadded. The resulting image (after cor- 
rection for exposure and smoothing with a a = 0.6' Gaussian) 
is shown in Fig.|| The elongation of the image and the varia- 
tion of the centroid of the X-ray emission on spatial scales of 
a few arcmin utes are consistent with those reported by A llen & 
Fabian (1994). To further emphasize the asymmetric structures 
in the surface brightness distribution, the image (/) was repre- 
sented as a sum of "symmetric" (/,) and "asymmetric" (I a ) com- 
ponents, where I s was calculated as the average surface bright- 
ness in 0.5' wide annuli centered at NGC4696. The asymmetric 
component is equal to the difference between the original im- 
age and the symmetric component (I a = I—I s ). Shown in Fig. || 
is the ratio of the asymmetric to the symmetric components, 
la/h = U -Is) /Us), smoothed with a a = 45" Gaussian. This 
image is, in fact, the relative deviation of the surface bright- 
ness from the average value at a given radius. Although such 
a decomposition of the surface brightness distribution into two 
components is of course somewhat arbitrary, there are several 



interesting features apparent in the figure. First, there is excess 
emission to the west of NGC4696, which is apparently respon- 
sible for the variation o f the centroid of the X- ray isophotes as 
found in the analysis of Allen & Fabian (1994). Note that in the 
raw image, the surface brightness peak coincides nicely with 
the optical position of NGC4696 (at least better than — 15" ). 
Second, there are two very clear excesses of emission. One 
lies close to the position of the SO galaxy NGC4696B (see also 
Allen & Fabian 1994) and another is centered on NGC4709 
(which was not detected by Allen & Fabian (1994)| who used 
cleaned data from only one pointing having a total exposure 
time of 3311 seconds). From Fig. [| it is clear that NGC4709 
marks the center of extended structure, lying to the east (and 
south-east) of NGC4696. Given the fact that NGC4709 is the 
dominant galaxy in the low velocity dispersion subgroup, Cen 
45 (see Lucey et al. 1986), this suggests that this extended emis- 
sion is associated with Cen 45. Of course, the intensity and ex- 
tent of the structure around NGC4709 depends on the assumed 
centroid position used for calculation of the "symmetric" com- 
ponent. However, the detection of the excess emission close to 
NGC4709 is stable against variations of the centroid position 
for the symmetric component. 

Since the Cen 45 group contains a factor of 2.5 fewer galaxies 
than Cen 30 (Lucey et al. 1986) and a significantly lower ve- 
locity dispersion (280 km s" 1 compared to 586 km s" 1 for the 
main group), a simple projection of the two subgroups would 
result in an excess surface brightness around NGC4709 (con- 
sistent with the ROSAT data), but a lower temperature in this 
region (contrary to the ASCA data, see Fig. [I]). We suggest that 
the most straight forward explanation would be to assume that 
the gas in the region around NGC4709 is heated due to the in- 
teraction between the two subgroups. Note that the direction 
of the line joining NGC4709 and NGC4696 is co-aligned with 
featur es at various spatial scales. As noted by A llen & Fabian 
(1994), the major axes of the elliptical X-ray isophotes around 
NGC4696 (at r> 100 arcsec) are aligned towards NGC4709 as 
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Fig. 3. — Left: The central (40 x 40 arcminutes) region of the ROSAT PS PC image corrected for vignetting overlayed onto optical 
image. The bright elliptical galaxy to the South East is NGC 4709; the galaxy to the North West is NGC 4696B. The X-ray image 
was smoothed with a a = 0.6' Gaussian. Right: Contour plot of the relative excess emission with respect to the azimuthally averaged 
surface brightness for the PSPC (i.e. Da '^ d ° e f d ). The image was smoothed with a a = 45" Gaussian. Contours start at 0.05 with a 
0.05 increment. The circle outlines the position of the PSPC support structure. Contours outside this circle were suppressed in order 
to avoid possible spurious features induced by the PSPC support structure. 



are the optical isophotes of NGC4696 itself. At scales up to a 
few degrees, this sa me orientation marks the chain of separate 
groups, identified by Lucey et al. (1986). Thus, it may be that 



the Centaurus cluster is accreting matter (in rather small incre- 
ments) along this preferential direction which perhaps marks 
the filamentary structure at large scales (see e.g. West, Jones, 
& Forman 1995). At the present time, the subgroup centered on 
NGC4709 is merging with the main cluster. 

Another interesting feature which can be seen in Fig. || is an 
X-ray source close to th e position of the SO galaxy NGC4696B. 



Allen & Fabian (1994) noted that this X-ray source probably is 



associated with NGC4696B, which is one of the three brightest 
galaxies in the field (only NGC4696 and NGC4709 are brighter 
than NGC4696B in the optical band). However, the X-ray po- 
sition is offset from the optical center of the galaxy by ~ 1' 
towards the north-east. Moreover, in Fig. [}]one can see a long 
and narrow filament which originates close to NGC4696B and 
extends to the north-east and then to the east. The filament is 
barely visible in the raw images and appears clearly only af- 
ter subtraction of the component centered on NGC4696. The 
significance of the excess at various positions in the feature (de- 
fined as the ratio of the signal at a given position to the statistical 
error associated with the smoothed image) corresponds roughly 
to 2-3 a. Note however that nearby pixels are not independent 
on the smoothed image. 

Of course it is difficult to prove that the observed excess is 
a single feature. One can speculate however on the the origin 
of this feature using the same line of a rguments as presented 
by Vikhlinin, Forman, & Jones (1996), who discussed a fila- 



ment of X-ray emission in Coma. The two simplest scenarios 
assume that such a feature may be due to cold gas (ram pres- 
sure stripped) from NGC4696B, like that observed for M86, 



or due to a tidally disrupted dark halo from a group of galax- 
ies during its passage through the central part of the Centaurus 
cluster. Future observations (AXAF, XMM and SPECTRUM- 
X-GAMMA) may distinguish between these two possibilities: 
a lower temperature and higher abundance of the excess emis- 
sion would suggest stripping, while a lack of spectral changes 
would argue in favor of a perturbation in the gravitational po- 
tential. 

The shift in the position of the peak of the X-ray emission 
from NGC4696B to the northeast of the optical position sup- 
ports the suggestion that the whole feature is somehow associ- 
ated with NGC4696B. Note also that the major axis of the op- 
tical emission of NGC4696B is also aligned in the direction of 
the tail. Assuming a radius for the feature of ^25 kpc, a length 
^200 kpc, and an enhanced electron density of ~ 10~ 3 cm" 3 
(required to produce the observed surface brightness increase), 
the total gas mass of the feature is of the order of 6 x 10 9 Mq. 
The required stripping rate can be estimated as 16 M Q V500 yr~ l , 
where V500 is the galaxy velocity in units of 500 km sec" 1 . This 
value is too high for steady mass ejection by this single galaxy, 
suggesting either that the situation is transient or a group of 
galaxies is responsible for this gas. Considering that Centau- 
rus is likely to be undergoing repeated mergers with smaller 
subgroups, we can expect that relatively gas rich galaxies are 
frequently entering the dense environment of this cluster. The 
evaporation time (assuming that this is cold gas stripped from 
NGC4696B) of such a feature due to thermal conduction would 
be short (~ 10 7 years) and one would have to assume that ther- 
mal conduction is suppressed below Spitzer's value, which is 
not unlikely (see e.g. Pistinner, L evinson, & Eichler 1996). As 
discussed by |Soker et al. (1991) , stripped cool gas from cluster 
galaxies may fall to the bottom of the potential well and provide 
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strong density fluctuations in the cooling flow region which are 
required to trigger thermal instabilities. It is interesting that 
the assumption that the observed filament arises from cool gas 
stripped from NGC4696B may indicate a rather high infall ve- 
locity for the cool gas. Indeed, if the bending of the tail towards 
NGC4696 is due to infall of the gas in the cluster potential, 
then the infall velocity should be comparable to the velocity 
of NGC4696B. The infall velocity of an overdense blob of ra- 
dius r cannot exceed its terminal velocity v, ~ vx y/(Apr)/(pR) 
(e.g. Nulsen 1986), where v# is the Keplerian velocity and R is 
the distance from the cluster center. If most of the mass in the 
feature is in X-ray emitting gas, then 



/Apr 
V PR 



1/2 



= V K 



( Ap 2 r 



-I 1 / 4 



(6) 



The term is the surface brightness enhancement towards 
the feature, which is roughly 0.2-0.3, while r/R ~ 0.1. There- 
fore, v, ~ 0.4v/f. Infall of the blob at this velocity may be 
consistent with the curvature of the feature, if the velocity of 
NGC4696B is not much higher than v^. Such a high infall ve- 
locity would result in the tail, arising from stripped gas (which 

is visible in the images, i.e. is not much less than unity), 
having one end lying near the cluster center. However, the 
tail may be disrupted as it moves through the cluster gas (e.g. 
Nulsen 1996), effectively decreasing the terminal velocity of 
the infalling gas as it fragments into smaller blobs. Considering 



typical parameters for cooling flow clusters, Soker et al. (1991) 
found an infall velocity for stable blobs of roughly 60 km s" 1 , 
assuming a magnetic field of ~ 1/iG. Such low infall velocities 
would mean that the shape of the tail is just a reflection of the 
projection of the galaxy trajectory (which therefore must have 
passed through the very center of the Centaurus cluster). 

4. DISCUSSION 

Despite of the large difference (~ 150 km s" 1 ) in the ve- 
locities of the Cen 30 and Cen 45 groups, ^ucey et al. (1986)| 



presented strong arguments that these two groups are located at 
the same distance. Luminosity functions, color-magnitude re- 
lations and galaxian-radius distributions support the assumption 
of a common distance and disagree with the contention that Cen 
45 is significantly more distant than Cen 30, in accordance with 
their line-of-sight velocities. Dressier (1995) came to the same 
conclusion studying surface brightness fluctuations of several 
galaxies from both groups. These results imply that the smaller 
Cen 45 group is merging with the dominant Cen 30 group. The 
analysis of X-ray data further supports this scenario. Indeed, if 
one were observing the projection of two independent clusters, 
we would expect the overlapping region to have a lower tem- 
perature due to the predicted lower temperature of the smaller 
group Cen 45. Instead, we find that the temperature is highest 
in the region close to the dominant galaxy of the Cen 45 group, 
which we interpret as heating due to the interaction of the two 
subgroups. Note however that the presence of the cooling flow 
around NGC4696 restricts our ability to correctly estimate the 
extent of the hot region. We might also expect that when a 
poor group is entering a dense cluster environment for the first 
time, the relatively cold gas will be stripped due to the high ram 
pressure. The marginally significant filament near NGC4696B, 
seen in the ROSAT PSPC images, might be explained as the 
result of such an interaction. 



Recently Stein, Jerjen, & Federspiel (1997) found, in the 
0.75°circle around NGC4696, a very large velocity dispersion 
for Cen 30 (933 ±118 km s" 1 ) and a very low dispersion for 
Cen 45 (131 ±43 km s" 1 ). The assumption that the main clus- 
ter has a velocity dispersion of 933 km s" 1 would make the so- 
called /3-problem very extreme. This parameter (which char- 
acterizes the ratio of energy per unit mass for galaxies to that 
in the gas) can be estimated from the X-ray surface brightness 
distribution, (3f, or calculated from the observed velocity dis- 
persion and gas temperature, (3 S . The slope of the PSPC surface 
brightness profile implies values of f3f ~ 0.4-0.45. However, 

P s = ^£f- ~1.4 for T ~ 4 keV and o> =933 km s _1 . Even 
though these values for (3 disagree for many clusters, for Cen- 
taurus this disagreement is far larger than usual. An alternative 
and more likely explanation for the large velocity dispersion is 
that a large scale filament, which defines the preferential direc- 
tion of accretion onto the cluster, lies close to the line of sight. 
This would explain the high observed velocity difference of the 
infalling group and wo uld also account for the large vel ocity 
dispersion observed by Stein, Jerjen, & Federspiel (1997) . The 
presence of a large scale filament crossing the center of Cen 30 
also would naturally explain the alignments already mentioned 
above. Thus, accretion along a filament, lying mostly along the 
line of sight, but also partially along the direction defined by 
the centers of Cen 30 and Cen 45, would explain the common 
alignments of 1) the major axis of NGC4696, 2) the filament 
emanating from NGC4696B and 3) the nearby chain of groups 
(Luceyetal. 1986). 

The effects of large scale filaments are seen around other 
clusters. For example, observations of A1689 suggest the pres- 
ence of a large scale filament directly along the line of sight 
(Daines et al. 1998). In A1689, the three brightest galaxies 
in the direction of the cluster lie within 20" of the cluster cen- 
ter, only one has a velocity equal to the cluster mean while the 
other two have velocity differences of +4767 km s" 1 and -2686 
km s -1 , far too large to be explained by the velocity dispersion 
of a relaxed ~ 10 keV temperature cluster. The nearby cluster 
A85 also shows the common alignment of several galaxy con- 
centrations/subclusters with the major axis of the central cD 
galaxy (Durret et al. 1997). The X-ray temperature map for 
A85 (Markevitch et al. 1998) shows a hot region, suggestive of 
a shock produced by one of the infalling subclusters. 

Large scale filaments and their effects on clusters also are 
found in recent numerical simulations which have adequate res- 
olution. Colberg et al. (1998) showed that clusters accrete 
matter from a few preferred directions, defined by filamentary 
structures, and that the accretion persists over cosmologically 
long times. Thus, it should not be surprising to see common 
alignments between central galaxy major axes, infalling sub- 
clusters, surrounding groups, and, even on the scale of neigh- 
boring clusters (West et al. 1995). With the availability of clus- 
ter temperature maps, the X-ray observations allow us to dif- 
ferentiate between superpositions (chance alignments) and true 
mergers where shock heating occurs. Hence, as shown in this 
paper for Centaurus, we can now study the details of the merger 
process which provide strong evidence for accretion along pre- 
ferred directions as predicted by numerical simulations of hier- 
archical structure formation. 
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